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Abstract 

We study the low-energy 3/2^ and l/2~ states of ^He and ^Li in a mi- 
croscopic cluster model. The scattering phase shifts of Bond (a -|- n) and 
of Schwandt (a + p), respectively, are reproduced well. We determine the 
resonance parameters by localizing the poles of the analytically continued 
S-matrix at complex energies. Our results differ from conventional i?-matrix 
resonance parameters, which were extracted from experimental data using the 
definition of a resonance based on the positions and widths of reaction cross 
section peaks. However, they nicely agree with the results of an extended 
i?-matrix method which works at complex energies. 

PACS numbers: 27.10.-hh, 24.30.Gd, 21.60.Gx 
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Wigner's i?-matrix is one of the most powerful tools in nuclear physics. It is rou- 
tinely used to analyze experimental data. Just to mention one example, the analysis of 
the ^^C(a, 7)^^0 process, which is a key reaction in nuclear astrophysics is heavily based 
on i?-matrix fits [@. The beauty of Wigner's method is that all scattering quantities are 
parametrized in terms of real, energy-independent quantities. It is true that the values of 
some of the parameters of the theory, namely the reduced- width amplitudes, 7ac, and energy 
eigenvalues. Ex, are somewhat arbitrary due to their dependence on the boundary- condition 
numbers, B^, and channel radii, Qc- However, any sufficiently robust description [i.e. one 
including enough levels) of the scattering data in a given energy region will give stable 
(independent of Oc, B^) resonance parameters when they are based on the actual complex- 
momentum poles and residues of the S'-matrix, as calculated from continuing the i?-matrix 
into the complex energy plane. 

We will call this method based on the actual S'-matrix pole structure in the complex 
plane the "extended" i?-matrix prescription for defining resonance parameters, in order to 
distinguish it from the usual resonance-parameter prescriptions that are defined entirely on 
the real energy axis. The real-energy parameters are easier to obtain, especially when they 
are extracted directly from the measurements, which of course exist only on the real energy 
axis of the physical sheet. However, it is sometimes difficult to interpret experimental results 
based solely on methods that work at real energies. We mention here two examples: the 
large cross section of the t{d,n)a reaction and the problem of the soft dipole resonance in 
neutron halo nuclei. Only analyses at complex energies were able to reveal that the large 
reaction cross section is caused by a shadow pole of the scattering matrix in the former case 
0,^, and that the soft dipole resonance does not exist in ^He in the latter one p. 

It is an intriguing question whether or not the results of methods that are not confined 
to real energies agree with those of conventional analyses for the relatively broad resonances 
encountered in light systems. As an exploratory investigation, in this paper we study the 
low-energy 3/2~ and l/2~ states of the ^He and ^Li nuclei. The parameters of these states 
have been determined from conventional i?-matrix analyses of certain experimental data. In 
those works the definition of the resonance position and width was based on the positions 
and widths of peaks in the reaction cross sections. 

For narrow isolated single channel resonances the two definitions of resonance parameters 
(S'-matrix pole and cross section peak, respectively) are consistent with each other, and 
give the same results. However, for broad resonances the results coming from the two 
definitions may be different, because only the scattering theoretical quantities defined at 
complex energies (e.g., the S'-matrix, the Fredholm determinant, and the Jost function on 
the multisheeted Riemann energy surface) contain the correct dynamical information. Even 
in the case of a narrow multichannel resonance the results coming from the two definitions 
can disagree, and studying the complex energy scattering quantities can give a much deeper 
insight into the dynamics of the problem 0,^. 

The extraction of the a + N resonance parameters from observables which are continued 
to complex energies was first suggested by Ahmed and Shanley . They pointed out that the 
determination of the resonance parameters from real energy observables is difficult because, 
e.g., the l/2~ phase shift does not even pass through 90 degrees. Here we determine these 
parameters from the analytic continuation of the a+N scattering matrix to complex energies 
in a microscopic model. We also extract the ^He and ^Li resonance parameters from the 
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extended i?-matrix method. 

Our model is a microscopic a + N Resonating Group Method (RGM) approach to the 
five-nucleon system. The trial function of the five-body problem reads 
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PaN, 



JM 



(1) 



where A is the intercluster antisymmetrizer, the PaN vector is the intercluster Jacobi co- 
ordinate, L and S is the total angular momentum and spin, respectively, and [...] denotes 
angular momentum coupling. While <l>^ (A^ = n or p) is a neutron or proton spin-isospin 
eigenstate, the antisymmetrized ground state {i = 1) and monopole excited states {i > 1) 
of the a particle are represented by the wave functions 



'ft' 



1,2, 
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Here (j)"^. is a translationally invariant shell-model wave function of the a particle with size 
parameter jSj and the Aij parameters are to be determined by minimizing the energy of 
the a particle Putting ([^) into the five-nucleon Schrodinger equation which contains a 
two-nucleon strong and Coulomb interaction, we arrive at an equation for the intercluster 
relative motion functions x- This equation is solved by utilizing a Kohn-Hulthen variational 
method for the S'-matrix, which uses square integrable basis functions matched with the 
correct scattering asymptotics p|. 

The input data for the i?-matrix studies are the cross sections and polarizations for all 
possible reactions involving ^He and ^Li. In the present RGM approach we concentrate on 
reproducing the a + N scattering phase shifts because they are most closely related to the 
usual definition of the resonance position and width. So, if our model reproduces the phase 
shifts, then our resonance parameters are hopefully close to the ones that characterize the 
reactions. Moreover, in order to analyze cross sections, we should build reaction mechanisms 
into the model, which would make this model rather phenomenological and ambiguous. 

We use the Minnesota effective N-N interaction 0. The same model (for the ^He and 
^Li subsystems) and interaction were used in [|l^,|ll|,^ to successfully describe the structure 
and beta delayed deuteron emission of ^He, and the three-body resonances of the A = 6 
nuclei. 

In Fig. 1 we show our S- and P-wave a + N phase shifts, together with the experimental 
data of Bond [O for a+n and Schwandt JTSl for a+p. A rather good agreement is observed, 



especially in the resonance region. We do not show the higher partial waves because they 
do not influence our results, and they are practically zero at low energies in agreement with 
the experiments. 

The experimental parameters of the low-lying ^He and ^Li states are listed in Table I. 
We compare our results with those of Barker and with the compilation |jl5|. The results 
of Ref. 1^ are also shown. 

The 3/2~ and l/2~ resonance parameters are determined by analytically continuing the 
S'-matrix to complex energies |^ . In practice this is done by solving the Schrodinger equation 
for the a — N relative motion at complex energies with the the following boundary condition 
for oo 
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Here and ki are the complex energies and wave numbers of the relative motions, and H~ 
and are the incoming and outgoing Coulomb functions, respectively. The function 5* 
has no physical meaning, except if it is singular at the energy e. Then S coincides with the 
physical S'-matrix describing a purely outgoing solution, that is a resonance. So we search for 
the poles of S at complex energies. All quantities are defined on the multisheeted Riemann 
energy surface, and are analytic almost everywhere. The complex Coulomb functions were 
calculated by using ||I6[| . 

The positions (Er) and widths (F) of the resonances are extracted from the e = E,. — 
iV /2 complex pole positions of the S'-matrix. We compare these parameters with those 
coming from the cross section peak definition of a resonance, used in conventional i?-matrix 
approaches. We can see in Table I, that our parameters are rather different from the R- 
matrix results of fl^. It is especially intriguing that the splitting of the 3/2 

M 



states is much smaller than in 



0. 



and 1/2- 

We explored the dependence of 



not to mention 

the phase shifts on the resonance parameters by slightly changing the N-N interaction. In 
Fig. 1(a) the dashed lines show the change of the resonant phase shifts when the resonance 
parameters of Table I are changed by 10%. We also checked the effect of including the d + t 
channel in Eq. (0), and found that the resonance parameters are little changed, provided 
the phase shifts have the same quality as in Fig. 1(a). 

We also used the extended i?-matrix method [§] to extract resonance parameters for the 
A = b ground- and first excited-states from multichannel i?-matrix analyses of reactions 
in the ^He and ^Li systems. The A = b analyses included the two-body channels N + a 
and d + t oi rf+^He, along with pseudo two-body configurations to represent the breakup 
channels n + p + toin + p+^He. Included in the n — a data are the differential elastic 
scattering cross sections of Morgan [0, Hoop [|I^], Niiler [^, and Shamu polarization 
and analyzing-power measurements by Sawers ||2l], Broste [22|, May [23|, and Perkins [2^|; 
and neutron total cross sections measured by Haesner The p — a data include the 

|2^, Garreta [|8[, and 



261, Jarmie 



differential elastic scattering cross sections of Freier 
Plattner |2^; polarization and analyzing-power measurements by Schwandt [|l^, Plattner 
and Hardekopf 



I0|; and polarization-transfer measurements by Keaton ||31|| . 
The results are given in Table I. It is remarkable that most of the resonance parameters 
nicely agree with the ROM results, but differ from the conventional i?-matrix results. Per- 
haps the only exception is the 1/2" state of ^Li, where the agreement between the ROM 
and extended i?-matrix results is not so good. The small differences between the extended 
i?-matrix and RGM results in the case of ^He probably come from the fact that our i?-matrix 



gives phase shifts that are slightly different from ||T3[ at higher energies. Our RGM results 
are also in good agreement with |^, where a higher order scattering amplitude expansion 
was used. 



In addition to the dependence of the resonance parameters of Ref. [0] on the mechanism 
by which ^He and ^Li are formed, there is also a marked dependence on channels radius. 
Barker has argued that this dependence can be used to determine a "best" value of channel 
radius, which for n -|- a is taken to be 5.5 fm. It is therefore quite interesting that when 
the n — a i?-function parameters from Table 8 of Ref. are used in the S'-matrix pole 



prescription P], the resulting resonance parameters, E-, 



3/2- 



0.77 MeV, F 



3/2- 



0.65 MeV; 
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El/2- = 2.10 MeV, ri/2- = 5.37 MeV, are in good agreement with the RGM values and 
with the extended i?-matrix values, defined for Qn-a = 3.0 fm. 

The simplest way to extract resonance parameters is to fit the cross section or phase 
shift data with Breit-Wigner forms. This is, however, an ambiguous procedure for broad 
resonances where the phase shift is not "ideal" {i.e., not going from to 180 degrees within 
a short energy interval). For an "ideal", isolated, narrow resonance the phase shift, given 
by scattering theory, behaves like tan 6{E) = 0.5r/{Er — E), which implies that d6/dE has 
a maximum a.t E = Er, and F = 2/ {d5 / dE)^^. This prescription to extract E,. and F is 
also used for broader resonances [^]. In Table II we show the 3/2~ and l/2~ parameters 
coming from this definition applied to the experimental and model phase shifts. One can 
see that this simple procedure provides resonance parameters which are close to the RGM 
results, except that it systematically overestimates the widths. 

In a very recent paper |3^ the authors extracted the resonance parameters of the ^He 
and ^Li ground states by determining the pole positions of the S'-matrix corresponding to 
3H(c/, 7)5He and ^}le{d, ^fU measurements |3|]. Their results, E3/2- (^He) = 0.8±0.02 MeV, 
F3/2-(^He) = 0.65 ±0.02 MeV, and E^/2-{^U) = 1.72 ±0.03 MeV, F3/2-(^Li) = 1.28 ±0.03 
MeV are in good agreement with our values in Table I. 

In summary, we have determined the parameters of the low-energy ^He and ^Li resonances 
from the complex pole positions of the a + N scattering matrix in a microscopic cluster 
model. Our results are different from the results coming from a conventional i?-matrix 
method, which define the resonance parameters based on the real energy properties of cross 
section peaks. However, they are in good agreement with the results of an extended i?-matrix 
method that works in the complex energy plane. We emphasize that the extended i?-matrix 
method involves no difference in the way that i?-matrix parameters are extracted from 
experimental data, but only in the way that they are subsequently used to define resonance 
parameters. However, it is quite clear that the dependence on the channel radius of the 
usual real-energy prescription for defining resonance parameters from i?-matrix parameters 
is characteristic of that prescription. 

The S'-matrix pole prescription gives consistent resonance parameters for the ground 
state and first excited state of the A = 5 nuclei, which are approximately independent of 
the method used to describe the nuclear dynamics or the reaction in which the resonance is 
observed, and in the case of the i?-matrix parametrizations, also independent of the channel 
radii and boundary conditions. We expect this would also be the case for relatively broad 
levels in other light systems, where different resonance-parameter prescriptions can lead to 
quite different results [^, and so we recommend using the complex 5'-pole prescription to 
specify resonance parameters in all cases. 
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FIGURES 

FIG. 1. Phase shifts for a + n (a) and a + p (b) scattering in our RGM model. The dashed 
lines are results from calculations where the resonance parameters differ from those in Table I by 
10%. Experimental data points are taken from (a) and (b), respectively. 20 degrees are 
added to the S1/2 phase shifts, for clarity. 
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TABLES 



TABLE L Parameters of the low-energy ^He and ^Li resonances in the center-of-mass frame, 
coming from the S'-matrix pole and cross section peak definitions, respectively. Er is the resonance 
position relative to the a + N threshold, and T is the full width at half maximum. All numbers 
are in MeV. 



^He 



Method 



Er[3/2-) r(3/2-) Er(l/2-) r(l/2^ Er[3/2~) r(3/2-) Er(l/2-) r(l/2-) 



Compilation 
_R-matrix, stripping 
_R-matrix, pickup ^4 f 
Scattering ampl. iGf 
5-matrix, RGM 
Extended ij-matrix 



0.89±0.05 

0.838±0.018 

0.869±0.003 

0.778 

0.76 

0.80 



0.60±0.02 

0.645±0.046 

0.723±0.019 

0.639 

0.63 

0.65 



4.89±1 

2.778±0.46 

3.449±0.4 

1.999 

1.89 

2.07 



4±1 
3.6±1.2 
5.3±2.3 
4.534 
5.20 
5.57 



1.96±0.05 

1.76±0.06 

1.86±0.01 

1.637 

1.67 

1.69 



Ri 1.5 
1.18±0.13 
1.44±0.08 
1.292 
1.33 
1.23 



7-12 
3.63±0.56 
4.54±0.5 
2.858 
2.70 
3.18 



5±2 
4.1±2.5 
6.1±2.8 
6.082 
6.25 
6.60 



TABLE II. Parameters of the low-energy ^He and ^Li resonances in the center-of-mass frame, 
determined by assuming that d5/dE has a maximum at E = Er, and F = 2 / [d5 / dE) ■ Er is the 
resonance position relative to the a + N threshold, and T is the full width at half maximum. All 
numbers are in MeV. 



'He 



'Li 



Phase shifts 



Er{2,/2^) r(3/2-) Er{l/2-) T{l/2-) Er{2,/2-) r(3/2-) Er{l/2-) T{l/2- 



Exp. data [|T2|,^ 

RGM 

i?-matrix 



0.77 


0.69 


2.13 


7.26 


1.53 


1.42 


2.77 


8.89 


0.76 


0.68 


2.07 


7.18 


1.67 


1.46 


2.92 


8.88 


0.75 


0.85 


2.21 


7.98 


1.67 


1.37 


3.35 


9.40 
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